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The global supersymmetry is spontaneously broken if and only if the ground-state energy is 
strictly positive. We propose to use this fact to observe the spontaneous supersymmetry breaking in 
euclidean lattice simulations. For lattice formulations that possess a manifest fermionic symmetry, 
there exists a natural choice of a hamiltonian operator that is consistent with a topological property 
of the Witten index. We confirm validity of our idea in models of the supersymmetric quantum 
mechanics. We then examine a possibility of a dynamical supersymmetry breaking in the two- 
dimensional J\f — (2, 2) super Yang- Mills theory with the gauge group SU(2), for which the Witten 
index is unknown. Differently from a recent conjectural claim, our numerical result tempts us to 
conclude that supersymmetry is not spontaneously broken in this system. 

PACS numbers: 11.15.Ha, ll.30.Pb, ll.30.Qc 



It is widely beUeved that supersymmetry is relevant 
in particle physics beyond the standard model and it is 
spontaneously broken by some mechanism. If supersym- 
metry is not spontaneously broken at the tree level of 
the loop expansion, it remains so to all orders of the 
loop expansion. There still exists, however, a possibil- 
ity that supersymmetry is spontaneously broken non- 
perturbatively. Precise study of such dynamical super- 
symmetry breaking remains elusive because we have no 
universal framework that defines supersymmetric (espe- 
cially gauge) theories at a non-perturbative level. 

Generally, the Witten index Tr(-l)^ Q, where F 
is the fermion number operator, provides an important 
clue. One can infer that the dynamical supersymmetry 
breaking does not occur in a wide class of supersymmetric 
models for which the Witten index can be computed to 
be non-zero. However, the Witten index is not a panacea. 
There exist physically interesting models for which it is 
very difficult to determine the Witten index and, in some 
cases, the index itself might be ill-defined due to a gap- 
less continuous spectrum. 

In this letter, we consider a possibility to observe the 
dynamical supersymmetry breaking in euclidean lattice 
simulations, in the light of recent developments on lattice 
formulation of supersymmetric theories 0, d, 0j S] ■ The 
conceptually clearest way to observe the spontaneous su- 
persymmetry breaking would be to examine the degen- 
eracy of bosons' and fermions' mass spectra in two-point 
correlation functions. Here, we propose an alternative 
method that is based on the following fact. The global 
supersymmetry is spontaneously broken if and only if 
the ground-state energy is strictly positive 6] . In princi- 
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pie, therefore, one can judge whether the supersymmetry 
breaking takes place or not if the ground-state energy can 
be computed. 

Let us first recall that the thermal average of the hamil- 
tonian H with the inverse temperature /3 is expressed by 
the euclidean functional integral as [2^ 

where S is the euclidean (i.e., imaginary-time) action 
and dfi symbolically denotes a functional integral mea- 
sure. In the right-hand side, the time-period of the 
system is taken to be /3. What is very important in 
Eq. ^ is the boundary condition in the temporal di- 
rection. It must be periodic for all bosonic variables and 
anti-periodic (aPBC) for all fermionic variables. (For all 
bosonic variables and for all variables with respect to 
spatial directions, we always assume the periodic bound- 
ary conditions.) In the large imaginary-time or the low- 
temperature limit f3 — + oo, only the ground state(s) con- 
tributes to Eq. H]) and the ground-state energy Eq is 
given by 

Eo = lim (i/)aPBC. (2) 

If Eq > 0, supersymmetry is spontaneously broken and 
it is not if Eq — 0. 

Suppose that in Eq. ((T|) one uses instead the periodic 
boundary condition (PBC) for all variables. Then the 
partition function is proportional to the Witten index 
A/pBc/pBcdA^e-^ = Tr(-l)^e-^« - Tr(-l)^ UM 
and Jp^^dfj, He^^ is proportional to the /3-derivative of 
the Witten index, that is always zero 

A/PBC / dfiHe'"^ = TY{-lfHe-''" = 0. (3) 

JPBC 

This independence of the Witten index on a parameter 
of the theory /3 is a consequence of the supersymmetry 
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algebra We also note that, when the Witten index 
is non-zero, Eq. ([3]) is not invariant under a shift of the 
origin of the energy H H + c. With the periodic 
boundary condition, wc thus have (independently of (3) 
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and clearly this provides no useful information on the 
ground-state energy. 

In Eq. ([1]), boundary conditions for realizing the ther- 
mal equilibrium explicitly breaks supersymmetry. In 
Eq. (121), we then observe how the effect of the tempera- 
ture (that is a conjugate variable to the energy) remains 
in the zero-temperature limit [3 ^ oo. If the effect re- 
mains, that is, i?o > in Eq. ([2]), we judge that the 
spontaneous supersymmetry breaking occurs. It is in- 
teresting to note that this procedure is quite analogous 
to a usual way to observe the spontaneous breaking of 
ordinary symmetries [1^. 

There are several issues to be clarified to embody the 
basic formula ([2]) in euclidean lattice formulation. First, 
the above argument assumes that a regularization to de- 
fine the functional integral does not break supersymme- 
try. Although lattice regularizations are generally irrec- 
oncilable with supersymmetry, for theories with extended 
supersymmetry, it is sometimes possible to set up a lat- 
tice regularization that preserves the invariance under a 
part of supersymmetry transformations [2, i, 4, 5] . Then, 
if the spacetime dimension is low enough, one may expect 
that the invariance under a full set of supersymmetry 
transformations is restored in the continuum limit. In 
what follows, we assume this sort of lattice regulariza- 
tion. 

Closely related to the above point, we have to properly 
choose a possible additive constant in the hamiltonian H . 
In other words, we have to correctly choose the origin of 
the energy. This is of course a crucially important point 
to judge the spontaneous supersymmetry breaking from 
the positivity of _Eo . A natural prescription to define the 
hamiltonian is to use the supersymmetry algebra [2^ . 
By the following reason, however, this issue of a "cor- 
rect" hamiltonian is somewhat delicate in the functional 
integral formulation based on the lagrangian. 

Suppose that the (for simplicity, off-shell) supersym- 
metry algebra is realized by the transformation law for 
variables appearing in the continuum lagrangian. This 
implies that there exists a fermionic transformation Q 
such that {Q,Q} = 2idQ, where Q is a fermionic trans- 
formation conjugate to Q and c^o is the time derivative. 
One would then expect from this algebra that the rela- 
tion iQQ ~ 2H holds, where Q is the Noether charge (in 
field theory, we use the Noether current instead) associ- 
ated with the transformation Q and H is the hamiltonian 
obtained from the lagrangian by the Legendre transfor- 
mation. 

In reality, this relation holds only up to equations of 
motion. Generally one ends up with iQQ/2 = H + 
(terms being proportional to equations of motion). The 



additional terms, that would be negligible in classical 
theory, cannot be neglected in general within the func- 
tional integral because those terms may give rise to con- 
tact terms at a coincident point, i.e., ultraviolet-divergent 
constants [25| . 

However, if a lattice formulation one adopts possesses 
at least one exactly-preserved fermionic symmetry, say 
the above Q, it is natural to adopt H = iQQ/2 as the 
definition of a hamiltonian. First, this structure is sug- 
gested from the supersymmetry algebra. Second, this 
choice has the correct origin of the energy in the sense 
that it is consistent with the topological property of the 
Witten index, Eq. ([3]), that is. 
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where we have assumed that the lattice action S and 
the integration measure d/i are invariant under the Q- 
transformation and the integral /pg^ d/i Q is finite. 
As already noted, when the Witten index is non-zero, this 
property fixes the origin of the energy uniquely. For these 
reasons, we adopt the Q-exactness of the hamiltonian, 
H = iQQ/2., as a working hypothesis in what follows. 

We examine our idea by applying it to a euclidean 
lattice formulation of the supersymmetric quantum me- 
chanics Q. The euclidean lattice action of this model 
can be taken as 0, [l^, Q| 



V(a;) {d + W"{4>{x)))^{x) 



W'{(f>{x))d(l){x] 



(6) 



where A = {x G aZ \ < x < (3} {a denotes the lattice 
spacing), d is the forward difference df{x) = f{x -I- a) — 
f{x) and V^a; = Z?) = = 0) for PBC and ipix = 

/?) — —%l]{x = 0) for aPBC. This lattice action is invariant 
under a lattice counterpart of one of the Af = 2 supersym- 
metry transformations, Q, that is defined by Q(j){x) = 
ijj{x), Qi}{x) = 0, Qi^{x) = F{x) - d<j){x) - W'{<j){x)) 
'AndQF{x) = d^{x)+W" {(l){x))^p{x). The corresponding 
continuum theory is invariant under also Q and the su- 
persymmetry algebra reads Q^ = Q =0 and {Q, Q} = 
—2d. Corresponding to this Q-symmetry, we have the 
Noether charge Q. By using a lattice transcription of this 
Noether charge Q{x) = —ipix) {id(j){x) — iW'{(l){x))) /a, 
we define the hamiltonian operator by H{x) = iQQ{x)/2. 
As elucidated above, this lattice hamiltonian has a cor- 
rect zero-point energy in the sense of Eq. ([5]). In the 
corresponding continuum theory Q , it is well-known that 
supersymmetry is spontaneously broken if and only if the 
number of zeros of the function W'{(f)) is even. 
As a definite example, we consider 

Wicl^ixj) = i(am)0(x)2 + ^{amf/' X<p{x)\ (7) 
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FIG. 1; lima^o(^^(2;))aPBc/?Ti as a function of the physical 
temporal size of the system /3m. The errors are only statistical 
ones. The dotted curve is the analytic expression for the A = 
case for which supersymmetry is not broken. 



FIG. 2: \ima-.o{H{x))s.PBc/'m and lima-,o{H{x))pBc/iri, as 
a function of the physical temporal size of the system /3m. 
The errors are only statistical ones. 



where m is a parameter which has the mass dimension 1. 
Supersymmetry is dynamically broken in the correspond- 
ing continuum theory when A 7^ 0. 

Fig. 1 is a result of Monte Carlo simulation for A = 
10. The continuum limit of the expectation value of 
the hamiltonian with the anti-periodic boundary condi- 
tion, limci_^o(^^(2;))aPBc/™i is plotted as a function of 
the physical temporal size of the system /3toJ2^. For 
p-m > 1, we have limQ^o(i?(a;))aPBc/™ — 1-1 ^27]. From 
this, we infer that supersymmetry is spontaneously bro- 
ken and this is indeed the right answer. We also numeri- 
cally observed that, in the present system, {H{x))p-BC is 
not well-defined, being consistent with Tr(— 1)-'^ = in 
the target theory (recall Eq. (g])) [TH - 

Next, as an example in which supersymmetry is not 
spontaneously broken, we consider 

WWx)) = l{amfcj,(x)\ (8) 

In this case, we numerically observed that both 
{H{x))pBC and (^f(x))aPBC are well-defined and, in 
Fig. 2, we plotted the continuum limit of these quan- 
tities as a function of /3m [H]- The figure shows that, 
in this case, lima^o(-ff(2;))pBC is consistent with zero for 
all temporal sizes (recall Eq. (|4]); in the target theory 
Tr(— 1)^ = 1) and limQ^o(^^(a;))aPBC approaches zero 
as the temporal size of the system is increased. From 
Eq. we conclude that Eq — within the error and 
supersymmetry is not broken. 

Having observed that our method works perfectly in 
the supersymmetric quantum mechanics, we now study 
the two-dimensional Af = (2, 2) super Yang-Mills theory 
by using a lattice formulation proposed in Ref. [l^. (See 
also Ref. [l7j.) For this seemingly simple supersymmet- 
ric system, the value of the Witten index and whether 
supersymmetry is spontaneously broken or not are not 



known, due to non-compact flat directions of the clas- 
sical potential energy. In fact, the authors of Ref. 
conjectured the dynamical supersymmetry breaking in 
this system with the gauge group SU(A'^c)- 

We numerically studied only the case of 
the gauge group SU(2). The physical size 

of our two-dimensional lattice A is (3 x L; 
A = {x G al? I < a;o < /3, Q<xi < L] and the 
boundary condition of a generic fermionic field ip is set 
to be, ■0(2;o = PtXi) = +^(a;o = 0,a;i) for PBC and 
V-Cxo = (3,xi) = -^{xq = Q,xi) for aPBC [29|. 

The point is that the lattice action and the inte- 
gration measure of the lattice formulation of Ref. [l^ 
are manifestly invariant under a lattice counterpart of 
a part of the M = (2, 2) supersymmetry transforma- 
tions, Q. By a similar reasoning as above, a hamilto- 
nian density T-L{x) is then defined by n{x) = QJ^{x)/2, 
where Jq{x) is a lattice transcription of the Noether 
current associated with another fermionic symmetry of 
the target continuum theory, Qo- This definition is con- 
sistent with the supersymmetry algebra in the twisted 
spinor basis, {Q,Qo} = 2ii9o. From the Q-invariance 
of the lattice action and of the integration measure, we 
have Jpgp d/iH(a;) e^"^ = (assuming that the integral 
J-p-QQ A^j.J§{x) e^^ is finite) that is analogous to Eq. (O. 
Then the ground-state (vacuum) energy density Sq is 
given by lim^^oo lima^o(W(a;))aPBC = ^o- We judge 
that the dynamical supersymmetry breaking takes place 
if £0 > and it does not ii Eq — 0. 

Our algorithm and the simulation code^that were 
developed by using FermiQCD/MDP [H, are al- 
most identical to those of Ref. [2l|. We use the hy- 
brid Monte Carlo algorithm to generate configurations 
in the quenched approximation. The effect of dynam- 
ical fermions is then afterward taken into account by 
reweighting configurations by the pfaffian of the Dirac 
operator (we do not introduce any mass terms of fermions 
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or bosons that would explicitly break the Q-symmetry). 
Although this is certainly a brute force method compared 
to a standard pseudo-fermion algorithm, its implementa- 
tion is much simpler and the validity has been observed 
for one-point Ward-Takahashi identities |21|]. 

The number of statistically- independent configurations 
we used is summarized in Table f , where Nt and Ns are 
the number of lattice points for the temporal and the 
spatial directions, respectively. The physical size of the 
spatial direction is fixed to be Lg = y/2. We used the 
cold start and set all scalar fields to be zero at the initial 
configuration. As initial thermalization, we discarded the 
first 10^ trajectories and then we stored configurations at 
each 10^ trajectories (the auto-correlation time was 10- 
20 trajectories). 

The result of our Monte Carlo simulation is Fig. 3. The 
result was obtained, as Fig. 4, by an extrapolation to the 
continuum a = by a linear x^-fit. In Fig. 3, we ob- 
serve that £q is consistent with zero within the error. We 
regard this as an indication that supersymmetry is not 
dynamically broken in the two-dimensional Af = (2, 2) 
super Yang-Mills theory with the gauge group SU(2). Of 
course, errors in our present result are still large and 
we cannot completely exclude a possibility of the super- 
symmetry breaking of 0(1) in £o/g^. Further reduction 
of statistical errors will allow us to conclude whether 
the scale of the dynamical supersymmetry breaking is 
0(1) or not [30]. If the above observation of unbro- 
ken supersymmetry is true, the expectation value of the 



hamiltonian density under the periodic boundary condi- 
tion {'H{x))pBc should be well-defined and vanishes for 
all (3. Actually, the real part of the expectation values 
Re(7i(a;))pBc for various /3g we plotted in Fig. 3 show 
that this requirement is met within errors. (The imagi- 
nary part Im(7i(a::))pBC is also consistent with zero and 
the errors are quite smaller than those for the real part.) 
This provides another support for our observation. 

To our knowledge, this work is the first instance that 
the dynamical supersymmetry breaking in a gauge field 
theory (for which the Witten index is not known) is inves- 
tigated by numerical simulation. It should be interesting 
to consider applications of the present method to other 
supersymmetric theories and gain an insight on a possi- 
ble supersymmetry breaking that is difficult to obtain in 
other ways. 
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TABLE I: The number of statistically-independent configurations we used. 
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am = 0.1, 0.05 and 0.02. We used 10* statistically inde- 
pendent configurations for each set of parameters. 
[27] In Fig. 1, we plotted also the exact ground-state energy 
Eo/m = 1.27616 that was obtained by a numerical di- 
agonalization of the corresponding hamiltonian. The dis- 
crepancy of our Monte Carlo result for l3m > 1 with this 
exact result can be understood by a systematic error as- 
sociated with the linear extrapolation to the continuum. 
[28] The result was obtained by an extrapolation to the con- 
tinuum a = by a linear of data computed at 
am — 0.1 and 0.05. We used 10* statistically-independent 
configurations for each set of parameters. 
[29] Space does not permit to reproduce relevant mathemati- 
cal expressions and full details of the simulation. We refer 
the reader to Ref. [l^ for these and for a further list of 
references. 

[30] For the present lattice model, we are at present develop- 
ing a simulation code with the pseudo-fermion and the 
RHMC algorithm. We hope this enable us to reduce the 
statistical errors, without increasing the number of con- 
figurations substantially, 
statistical ones. 
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FIG. 4: Linear x^-extrapolation of Re(7i(a;))aPBc/<7^ to the 
continuum a = for various values of (3g. The errors are only 
statistical ones. 



